addpath(genpath('BlochSimulationCode'))
%% 
clear all
alpha1= [1:50];
clear NSRF10
for jj = 1:length(alpha1)
    T1_all = 1:1:1700;M0 = 1;
    TR=3.15*10;alpha=alpha1(jj)/180*pi;M0=1;
    NSRF10 = M0*(1-exp(-TR./T1_all))./(1-cos(alpha).*exp(-TR./T1_all));
    NSRF10_all(:,jj) = NSRF10;
    
    TR=3.15*7;
    NSRF7 = M0*(1-exp(-TR./T1_all))./(1-cos(alpha).*exp(-TR./T1_all));
    NSRF7_all(:,jj) = NSRF7;

    TR=3.15*4;
    NSRF4 = M0*(1-exp(-TR./T1_all))./(1-cos(alpha).*exp(-TR./T1_all));
    NSRF4_all(:,jj) = NSRF4;
end


figure; imagesc(abs(NSRF10_all),[0 1]); axis square;
set(gca,'fontweight','bold','FontSize',25,'LineWidth',3)
set(gca, 'tickLength',[0;0])

figure; imagesc(abs(NSRF7_all),[0 1]); axis square;
set(gca,'fontweight','bold','FontSize',25,'LineWidth',3)
set(gca, 'tickLength',[0;0])

figure; imagesc(abs(NSRF4_all),[0 1]); axis square;
set(gca,'fontweight','bold','FontSize',25,'LineWidth',3)
set(gca, 'tickLength',[0;0])

figure; imagesc((NSRF10_all-NSRF4_all),[0 .2]); axis square;
set(gca,'fontweight','bold','FontSize',25,'LineWidth',3)
set(gca, 'tickLength',[0;0])

figure; imagesc((NSRF7_all-NSRF4_all),[0 .2]); axis square;
set(gca,'fontweight','bold','FontSize',25,'LineWidth',3)
set(gca, 'tickLength',[0;0])




